! ******************************
FINISH
/CLEAR
/UNITS,SI 
/FILN,E101
/TITLE, MECHANICAL STRESS ANALYSIS OF FIXED TUBE-SHEET EXCHANGER
! ***************β趨***************
*SET,TSHELL,17          !ǳͲں
*SET,HSHELL,2500        
!ǳͲһԣؿܰʱԳԽģҿǳ̷ӦΪǳܳȵһ롣
*SET,DSHELL,806         !ǳͲھ
*SET,HPSH,20            !ǳ̲ܰ͹߶
*SET,RSHELL,15          !ǳ̲ܰԲǰ뾶
*SET,SHIGH,100          !ܰ
*SET,TCHANNEL,18.5      !ܳͲں
*SET,HCHANNEL,350       !ܳͲ
*SET,HPCH,20            !̲ܳܰ͹߶
*SET,RCHANNEL,15        !̲ܳܰԲǰ뾶
*SET,DP,25              !ȹֱ
*SET,SP,38              !ȹܹܼ
*SET,TPIPE,2.5          !ȹܱں
*SET,TINSULAT,50        !²
*SET,NY,20              !ȹ
*DIM,NXZ,,NY            !һάÿдŸлȹܸ
*SET,NXZ(1),6,6,6,6,5,6,5,6,5,5
*SET,NXZ(11),5,5,4,5,4,4,3,3,2,1
! ***************ʷֲ趨***************
*SET,NGH,8              !ȹʷ
*SET,NGMESH,40          !ȹʷ
*SET,DENSITYMESHG,2     !ȹʷּܱ
*SET,NR0,60             !ܰ廷ʷ
*SET,NHOU,2             !ǳ̺͹ܳͲںʷ
*SET,NINSULAT,3         !²ںʷ
*SET,NHSHELL,40         !ǳͲʷ
*SET,NCHANNEL,8         !ܳͲʷ
*SET,NR,4               !ܰԲʷ
*SET,NSHIGH,4           !ܰʷ
*SET,NHPSH,2            !ǳ̲ܰ͹߶ʷ
*SET,NHPCH,2            !̲ܳܰ͹߶ʷ
! ***************غ***************
*SET,PS,0.58            !ǳѹ
*SET,PT,2               !ܳѹ
! ***************ǰ***************
/PREP7
ET,1,SHELL63            !帨ģĿǵԪ
ET,2,SOLID45            !Ӧ嵥Ԫ
MP,EX,1,2.1E5           !ϵĵģ,1Ų϶ӦͲְ
MP,NUXY,1,0.3           !ϵĲɱ
MP,EX,2,202780          !ϵĵģ,2Ų϶Ӧְ
MP,NUXY,2,0.3           !ϵĲɱ
MP,EX,3,2E5             !ϵĵģ,3Ų϶Ӧܰͼ
MP,NUXY,3,0.3           !ϵĲɱ
MP,EX,4,93710           !ϵĵģ,4Ų϶Ӧȹֹܸ
MP,NUXY,4,0.3           !ϵĲɱ
! ***************ģ***************
!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!ܰ±,ķ֮һ
*IF,DSHELL/2.0D0-RSHELL,LE,DSHELL/2.0D0+TSHELL-TCHANNEL-RCHANNEL,THEN !жϲȷһܰԲƽ뾶ԲǱԵСΪ5MM
ROL=DSHELL/2.0D0-RSHELL-5
*ELSE
ROL=DSHELL/2.0D0+TSHELL-TCHANNEL-RCHANNEL-5
*ENDIF
PCIRC,,ROL,,90          !ɹܰԲƽ棬Ϊ1
*AFUN,DEG               !תɽǶ
NTOL=0                  !ȹ
*DO,J,1,NY              !ѭԻȹֱΪֱԲƽ
NTEMP=J
*DO,I,1,NXZ(NTEMP)
*IF,NTEMP,NE,NINT(NTEMP/2)*2,THEN                     !ݵǰȹżȷƫƵԲ
WPOFFS,SP*SIN(60)+2*SP*SIN(60)*(I-1),SP*(J-1)*COS(60) !JΪʱеIȹԲ
*ELSE
WPOFFS,2*SP*SIN(60)*(I-1),SP*(J-1)*COS(60)            !JΪżʱеIȹԲ
*ENDIF
PCIRC,,DP/2.0D0         !ԻȹֱΪֱԲƽ
NTOL=NTOL+1             !1
WPCSYS,1,0              !Բ
*ENDDO
*ENDDO
ASEL,U,AREA,,1          !ѡгܰԲƽ
CM,BASE1,AREA           !һѡеԻȹֱΪֱԲƽΪBASE1
ALLSEL                  !ȫѡ
ASBA,1,BASE1            !㣬ܰ±
!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!ȹ±,ķ֮һ
*DO,J,1,NY              !ѭȹ±Բ棬һѭͬ
NTEMPJ=J
*DO,I,1,NXZ(NTEMPJ)
NTEMPI=I
*IF,NTEMPJ,NE,NINT(NTEMPJ/2)*2,THEN
WPOFFS,SP*SIN(60)+2*SP*SIN(60)*(I-1),SP*(J-1)*COS(60)
*ELSE
WPOFFS,2*SP*SIN(60)*(I-1),SP*(J-1)*COS(60)
*ENDIF
*IF,NTEMPJ,EQ,1,THEN
PCIRC,DP/2.0D0-TPIPE,DP/2.0D0,,90
PCIRC,DP/2.0D0-TPIPE,DP/2.0D0,90,180
*ELSE
*IF,NTEMPJ,EQ,NINT(NTEMPJ/2)*2,THEN
*IF,NTEMPI,EQ,1,THEN
PCIRC,DP/2.0D0-TPIPE,DP/2.0D0,270,360
PCIRC,DP/2.0D0-TPIPE,DP/2.0D0,,90
*ELSE
PCIRC,DP/2.0D0-TPIPE,DP/2.0D0 
*ENDIF
*ELSE
PCIRC,DP/2.0D0-TPIPE,DP/2.0D0 
*ENDIF
*ENDIF
WPCSYS,1,0
*ENDDO
*ENDDO
AGLUE,ALL               !²ĻȹԲܰƽճΪһ
NUMCMP,AREA             !ʹŴС
NUMCMP,LINE             !ʹߵߺŴС
LSEL,S,RADIUS,,ROL      !ѡйܰԲƽԵ
ASLL,S,0                !ѡص棬ܰ±
CM,BASPLT,AREA          !ܰ±ΪBASPLT
ALLSEL                  !ȫѡ
ASEL,U,AREA,,BASPLT     !ѡлȹ±Բ
CM,BASTUB,AREA          !ȹ±ԲΪBASTUB
ALLSEL                  !ȫѡ
!!!!!!!!!!!!!!!!!!!ʷֲ
LSEL,S,RADIUS,,ROL      !ѡйܰԲƽԵ
LESIZE,ALL,,,NR0        !ָߵʷ
ALLSEL                  !ȫѡ
LSEL,S,RADIUS,,DP/2.0D0-TPIPE,DP/2.0D0  !ѡйɻȹ±ԲԲ
LESIZE,ALL,,,NGH/4      !ָЩߵʷ
ALLSEL                  !ȫѡ
LSEL,S,LENGTH,,TPIPE    !ѡгΪȹܱں
LESIZE,ALL,,,1          !ȹܱںʷ1
ALLSEL                  !ȫѡ
!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!ܰ
MSHAPE,0,2D             !ָ浥Ԫ״Ϊı
TYPE,1                  !ָ浥Ԫ
ESIZE,SP-DP             !嵥Ԫߴ
MAT,3                   !趨
AMESH,BASPLT            !Բܰ±񻮷
K,10001,0,0,0           !첼ܰ±Ĺؼ
K,10002,0,0,TCHANNEL    !첼ܰ±Ĺؼ
K,10003,0,0,SHIGH-TSHELL!첼ܰ±Ĺؼ
K,10004,0,0,SHIGH       !첼ܰ±Ĺؼ
NUMCMP,LINE             !ʹߵߺŴС
L,10001,10002,NHOU      !ɲܰ±츨ߣָʷ
L,10002,10003,NSHIGH    !ɲܰ±츨ߣָʷ
L,10003,10004,NHOU      !ɲܰ±츨ߣָʷ
*GET,LINE0,LINE,,NUM,MAX!õߵߺ
TYPE,2                  !ָ嵥Ԫ
VDRAG,BASPLT,,,,,,LINE0-2,LINE0-1,LINE0 !ܰ嵥Ԫ
ACLEAR,BASPLT           !ܰ±浥Ԫ
!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!ȹ
MSHAPE,0,2D             !ָ浥Ԫ״Ϊı
TYPE,1                  !ָ浥Ԫ
MAT,4                   !趨
AMESH,BASTUB            !Իȹ±Բ񻮷
K,10005,0,0,SHIGH+RSHELL+HPSH+HSHELL !컻ȹ±ԲĹؼ
L,10004,10005           !ɻȹ±Բ츨
*GET,LINE1,LINE,,NUM,MAX!õߵߺ
LESIZE,LINE1,,,NGMESH,DENSITYMESHG  !ָʷַͼܱ
TYPE,2                  !ָ嵥Ԫ
VDRAG,BASTUB,,,,,,LINE0-2,LINE0-1,LINE0,LINE1 !ȹܵ嵥Ԫ
ACLEAR,BASTUB           !ȹ±Բ浥Ԫ
!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!ܳͲ
WPROTA,,90              !ת
NUMCMP,AREA             !ʹŴС
*GET,NAREA1,AREA,,NUM,MAX                    !õǰ
RECTANG,DSHELL/2.0D0+TSHELL-TCHANNEL,DSHELL/2.0D0+TSHELL,-(RCHANNEL+HPCH+HCHANNEL),-(RCHANNEL+HPCH) !ɹܳͲ
LSEL,S,LOC,Z,-(RCHANNEL+HPCH+HCHANNEL/2.0D0) !ѡ߶
LESIZE,ALL,,,NCHANNEL                        !ָ߶ʷ
ALLSEL                                       !ȫѡ
!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!ܰǳͲ͹ܳͲ
RECTANG,DSHELL/2.0D0+TSHELL-TCHANNEL,DSHELL/2.0D0+TSHELL,-(RCHANNEL+HPCH),-RCHANNEL  !ɹ̲ܰܳ͹
LSEL,S,LOC,Z,-(RCHANNEL+HPCH/2.0D0)          !ѡ߶
LESIZE,ALL,,,NHPCH                           !ָ߶ʷ
ALLSEL                                       !ȫѡ
NUMCMP,LINE                                                !ʹߵߺŴС
*GET,LINE2,LINE,,NUM,MAX                                   !õǰߺ
K,200001,DSHELL/2.0D0+TSHELL-TCHANNEL-RCHANNEL,0,-RCHANNEL !̲ܰܳԲԲ
K,400001,DSHELL/2.0D0-RSHELL,,SHIGH+RSHELL
K,300001,DSHELL/2.0D0+TSHELL-TCHANNEL,0,-RCHANNEL          !̲ܰܳԲǻ
K,300002,DSHELL/2.0D0+TSHELL-TCHANNEL-RCHANNEL+RCHANNEL*COS(45),,-RCHANNEL+RCHANNEL*COS(45)!̲ܰܳԲǻе
K,300003,DSHELL/2.0D0+TSHELL-TCHANNEL-RCHANNEL,,,          !̲ܰܳԲǻյ
K,300004,ROL,,,                              !γɴĹؼ
K,300005,ROL,,TCHANNEL                       !γɴĹؼ
K,300006,ROL,,SHIGH-TSHELL                   !γɴĹؼ
K,300007,ROL,,SHIGH                          !γɴĹؼ
K,300008,DSHELL/2.0D0-RSHELL,,SHIGH          !γɴĹؼ
K,300009,DSHELL/2.0D0-RSHELL+RSHELL*COS(45),,SHIGH+(RSHELL-RSHELL*COS(45)) !γɴĹؼ
K,300010,DSHELL/2.0D0,,SHIGH+RSHELL          !γɴĹؼ
K,300011,DSHELL/2.0D0+TSHELL,,SHIGH+RSHELL   !γɴĹؼ
K,300012,DSHELL/2.0D0+TSHELL,,SHIGH-TSHELL   !γɴĹؼ
K,300013,DSHELL/2.0D0+TSHELL,,TCHANNEL       !γɴĹؼ
K,300014,DSHELL/2.0D0+TSHELL,,-RCHANNEL      !γɴĹؼ
K,300015,DSHELL/2.0D0+TSHELL-TCHANNEL-RCHANNEL,,TCHANNEL !γɴĹؼ
K,300016,DSHELL/2.0D0-TSHELL,,SHIGH-TSHELL   !γɴĹؼ
LARC,300001,300002,200001,RCHANNEL           !ɴԵLINE2+1
LESIZE,LINE2+1,,,NR/2                        !ָʷ
L,300001,300014,NHOU                         !񻮷ָLINE2+2
L,300002,300013,NHOU                         !񻮷ָLINE2+3
L,300013,300014,NR/2                         !ɴԵLINE2+4
AL,LINE2+1,LINE2+2,LINE2+3,LINE2+4           !ɴ
LARC,300002,300003,200001,RCHANNEL           !ɴԵLINE2+5
LESIZE,LINE2+5,,,NR/2                        !ָʷ
L,300003,300015,NHOU                         !񻮷ָLINE2+6
L,300013,300015,NR/2                         !񻮷ָLINE2+7
AL,LINE2+3,LINE2+5,LINE2+6,LINE2+7           !ɴ
L,300003,300004,1                            !ɴԵLINE2+8
L,300004,300005,NHOU                         !񻮷ָLINE2+9
L,300005,300006,NSHIGH                       !񻮷ָLINE2+10
L,300006,300007,NHOU                         !񻮷ָLINE2+11
L,300007,300008,1                            !ɴԵLINE2+12
L,300008,300016,NHOU                         !񻮷ָLINE2+13
L,300016,300015,NSHIGH                       !񻮷ָLINE2+14
AL,LINE2+6,LINE2+8,LINE2+9,LINE2+10,LINE2+11,LINE2+12,LINE2+13,LINE2+14 !ɴ
L,300012,300016,NR/2                         !񻮷ָLINE2+15
L,300012,300013,NSHIGH                       !ɴԵLINE2+16
AL,LINE2+7,LINE2+14,LINE2+15,LINE2+16        !ɴ
L,300009,300012,NHOU                         !񻮷ָLINE2+17
LARC,300008,300009,400001,RSHELL             !ɴԵLINE2+18
LESIZE,LINE2+18,,,NR/2                       !ָʷ
AL,LINE2+13,LINE2+15,LINE2+17,LINE2+18       !ɴ
LARC,300009,300010,400001,RSHELL             !ɴԵLINE2+19
LESIZE,LINE2+19,,,NR/2                       !ָʷ
L,300010,300011,NHOU                         !񻮷ָLINE2+20
L,300011,300012,NR/2                         !ɴԵLINE2+21
AL,LINE2+17,LINE2+19,LINE2+20,LINE2+21       !ɴ
RECTANG,DSHELL/2.0D0,DSHELL/2.0D0+TSHELL,SHIGH+RSHELL,SHIGH+RSHELL+HPSH !ɹ̲ܰܳ͹
LSEL,S,LOC,Z,SHIGH+RSHELL+HPCH/2.0D0         !ѡ߶
LESIZE,ALL,,,NHPSH                           !ָ߶ʷ
ALLSEL                                       !ȫѡ
!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!ǳͲ
RECTANG,DSHELL/2.0D0,DSHELL/2.0D0+TSHELL,SHIGH+RSHELL+HPSH,SHIGH+RSHELL+HPSH+HSHELL
LSEL,S,LOC,Z,SHIGH+RSHELL+HPCH+HSHELL/2.0D0     !ѡ߶
LESIZE,ALL,,,NHSHELL                            !ָ߶ʷ
ALLSEL                                          !ȫѡ
LSEL,S,LOC,X,DSHELL/2.0D0+TSHELL-TCHANNEL/2.0D0 !ѡ߶
LESIZE,ALL,,,NHOU                            !ָ߶ʷ
ALLSEL                                       !ȫѡ
LSEL,S,LOC,X,DSHELL/2.0D0+TSHELL/2.0D0       !ѡ߶
LESIZE,ALL,,,NHOU                            !ָ߶ʷ
ALLSEL                                       !ȫѡ
!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
MSHAPE,0,2D                                  !ָ浥Ԫ״Ϊı
TYPE,1                                       !ָ浥Ԫ
MAT,2                                        !趨
AMESH,NAREA1+1                               !Թ񻮷
TYPE,2                                       !ָ嵥Ԫ
ESIZE,,NR0                                   !嵥Ԫߴ
VROTAT,NAREA1+1,,,,,,10001,10002,90,1        !嵥Ԫ
!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!ܰ
MSHAPE,2,2D                                  !ָ浥Ԫ״Ϊı
TYPE,1                                       !ָ浥Ԫ
MAT,3                                        !趨
AMESH,NAREA1+2,NAREA1+9                      !Թܰ񻮷
ESIZE,,NR0                                   !嵥Ԫߴ
VROTAT,NAREA1+2,NAREA1+3,NAREA1+4,NAREA1+5,,,10001,10002,90,1 !ܰ嵥Ԫ
VROTAT,NAREA1+6,NAREA1+7,NAREA1+8,NAREA1+9,,,10001,10002,90,1 !ܰ嵥Ԫ
!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!Ͳ
MSHAPE,2,2D                                  !ָ浥Ԫ״Ϊı
TYPE,1                                       !ָ浥Ԫ
MAT,1                                        !趨
AMESH,NAREA1+10                              !Ͳ񻮷
ESIZE,,NR0                                   !嵥Ԫߴ
VROTAT,NAREA1+10,,,,,,10001,10002,90,1       !Ͳ嵥Ԫ
ACLEAR,NAREA1+1,NAREA1+10,1                  !浥Ԫ
NUMMRG,ALL                                   !غϵĵ㡢ߡϢϲ
NUMCMP,ALL                                   !е㡢ߡϢ
FINISH
! ******************************
/SOLU
NSEL,S,LOC,Z,SHIGH+RSHELL+HPSH+HSHELL        !ѡ
D,ALL,UZ,0                                   !Լ
ALLSEL                                       !ȫѡ
NSEL,S,LOC,X,0                               !ѡ
D,ALL,UX,0                                   !Լ
ALLSEL                                       !ȫѡ
NSEL,S,LOC,Y,0                               !ѡ
D,ALL,UY,0                                   !Լ
ALLSEL                                       !ȫѡ
LSEL,S,RADIUS,,DSHELL/2.0D0                  !ѡ
LSEL,A,RADIUS,,DSHELL/2.0D0-RSHELL           !ѡ
LSEL,A,RADIUS,,ROL                           !ѡ
LSEL,U,LOC,Z,-RCHANNEL-HPCH-HCHANNEL,SHIGH-TSHELL !ɾ
ASLL,S,0                                     !ѡص
ASEL,U,LOC,Z,SHIGH-TSHELL/2.0D0              !ɾ
ASEL,U,LOC,Z,SHIGH+RSHELL                    !ɾ
ASEL,U,LOC,Z,SHIGH+RSHELL+HPSH               !ɾ
ASEL,U,LOC,Z,SHIGH+RSHELL+HPSH+HSHELL        !ɾ
SFA,ALL,,PRES,PS                             !Թܰϱʩѹغ
ALLSEL                                       !ȫѡ
LSEL,S,RADIUS,,DP/2.0D0                      !ѡ
LSEL,R,LOC,Z,SHIGH+RSHELL+HPSH+HSHELL        !ѡ
ASLL,S,0                                     !ѡص
ASEL,U,LOC,Z,SHIGH+RSHELL+HPSH+HSHELL        !ɾ
SFA,ALL,,PRES,PS                             !Իȹʩѹغ
ALLSEL                                       !ȫѡ
LSEL,S,RADIUS,,DP/2.0D0-TPIPE                !ѡ
ASLL,S,0                                     !ѡص
LSEL,ALL                                     !ѡ߶
LSEL,S,RADIUS,,DP/2.0D0                      !ѡ
ASLL,U,0                                     !ѡص
SFA,ALL,,PRES,PT                             !Իȹڱʩѹغ
ALLSEL                                       !ȫѡ
LSEL,S,RADIUS,,DSHELL/2.0D0+TSHELL-TCHANNEL  !ѡ
LSEL,A,RADIUS,,DSHELL/2.0D0+TSHELL-TCHANNEL-RCHANNEL !ѡ
LSEL,A,RADIUS,,ROL                           !ѡ
LSEL,U,LOC,Z,TCHANNEL,SHIGH+RSHELL+HPSH+HSHELL !ɾ
ASLL,S,0                                     !ѡص
ASEL,U,LOC,Z,TCHANNEL/2.0D0                  !ɾ
ASEL,U,LOC,Z,-RCHANNEL                       !ɾ
ASEL,U,LOC,Z,-RCHANNEL-HPCH                  !ɾ
ASEL,U,LOC,Z,-RCHANNEL-HPCH-HCHANNEL         !ɾ
SFA,ALL,,PRES,PT                             !Թܰ±͹ڱʩѹغ
ALLSEL                                       !ȫѡ
LSEL,S,RADIUS,,DP/2.0D0-TPIPE                !ѡ
LSEL,R,LOC,Z,0                               !ѡ
ASLL,S,0                                     !ѡص
ASEL,U,LOC,Z,TCHANNEL/2.0D0                  !ɾ
SFA,ALL,,PRES,PT                             !Իȹ±ʩѹغ
ALLSEL                                       !ȫѡ
*SET,PN,(DSHELL/2.0D0+TSHELL-TCHANNEL)**2*PT/((DSHELL/2.0D0+TSHELL)**2-(DSHELL/2.0+TSHELL-TCHANNEL)**2)   !ƽ
ASEL,S,LOC,Z,-RCHANNEL-HPCH-HCHANNEL         !ѡ
SFA,ALL,,PRES,-PN                            !ʩƽ
ALLSEL
SOLVE                                        !
FINISH
! ******************************
/POST1
PLNSOL,S,INT,0,1     	                     !Ӧǿͼ
! ***************Ӧ·1--1***************
*SET,FIRPOT,27934                            !ڱڽڵ
*SET,SECPOT,27932                            !ڽڵ
PATH,L1,2                                    !·Ƽ·
PPATH,1,FIRPOT                               !·
PPATH,2,SECPOT                               !·յ
PRSECT, ,0                                   !г·Ի
! ***************Ӧ·2--2***************
*SET,FIRPOT,27972                            !ڱڽڵ
*SET,SECPOT,27968                            !ڽڵ
PATH,L2,2                                    !·Ƽ·
PPATH,1,FIRPOT                               !·
PPATH,2,SECPOT                               !·յ
PRSECT, ,0                                   !г·Ի
! ***************Ӧ·3--3***************
*SET,FIRPOT,121                              !ڱڽڵ
*SET,SECPOT,7471                             !ڽڵ
PATH,L3,2                                    !·Ƽ·
PPATH,1,FIRPOT                               !·
PPATH,2,SECPOT                               !·յ
PRSECT, ,0                                   !г·Ի
FINISH
